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Abstract 

A calculation of the spectroscopic properties, energy levels and elec- 
r*j , tromagnetic transitions and moments, of the ten nuclei 128 > 130 Te, 128 > 130 i ; 

12S ' 130 Xe, 129 > 131 I, 12 7>i29 Te w ithin the framework of the interacting boson 
model (IBM-2) and its extensions (IBFM-2 and IBFFM-2) is presented. 
The wave functions so obtained are used to calculate single-/? and 2^/3/? 
Cn 1 matrix elements for 128 . 130 Te -> 128 ' 130 I _> 128 - 130 Xe decay. Use of the ef- 

fective value of the axial vector coupling constant gA,eS,p extracted from 
single-/? produces results for 2i/j3/3 in agreement with experiment. 



O '. 1 Introduction 

In recent years, the possible measurement of the absolute mass scale of neutrinos 
through 0^/3/3 decay has become of considerable interest. This decay takes place 
only when the neutrino is a Majorana particle with finite mass. Its occurrence 
has not been confirmed yet and is at the present time the subject of many 
experimental investigations. Concomitant with the Oj^/3/3 there is the 2vf3f3 
decay mode. This mode is allowed by the standard model and it has now been 
\ observed in several nuclei pQ. While the 0vj3/3 mode can be safely calculated 

in the closure approximation, since the average virtual neutrino momentum is 
| of order 100 MeV/c and thus well above the scale of nuclear excitations, the 

£f~) ■ closure approximation is not expected to be good for 2j//3/3 where the neutrino 

momentum is of order of few MeV/c and thus of the same scale of nuclear 
excitations. 

2v(3/3 decay without the closure approximation has been calculated within 
r> • the framework of QRPA [3] and LSSM |3j. In this paper, we initiate a new ap- 

\ proach to calculate 2vj3j3 without the closure approximation within the frame- 

work of the interacting boson model (IBM-2) and its extensions (IBFM-2 and 
IBFFM-2) [4]. This latter model has been used extensively to calculate spec- 
tra of odd-even medium mass and heavy mass nuclei (IBFM-2) and of odd-odd 
nuclei (IBFFM-2) [5], which are crucial for the calculation of 2^/3/3 decay. Af- 
ter a description of the IBFM formalism, we proceed to do a calculation of 
two-neutrino double-/? decays for 128 Te-> 128 Xe and 130 Te-> 130 Xe. The aim of 
the paper is two-fold: (i) first and foremost we want to understand what is 
the mechanism of 2^/3/3 in Te, that is, what intermediate states in the odd-odd 



1 This is a pre-copy-editing, author-produced PDF of an article accepted for publication 
in Progress of Theoretical and Experimental Physics following peer review. The definitive 
publisher-authenticated version will be available online. 



1 



nucleus contribute to the decay and (ii) from a comparison of our calculated 
matrix elements with experimental single-/3 and 2vj3j3 decay, extract the value 
of the effective axial vector coupling constant, gA,cS- 



2 Two-neutrino double-/? decay in IBFM 
2.1 Gamow- Teller and Fermi transitions 

The inverse half-life r _1 for 2v double- j3 decay has been derived by several 
authors [HI H E] - We use here the formulation of Tomoda [pj as adapted in [TO] . 
For transitions Of — > Oj, r _1 can be factorized to a good approximation as 

[rffM -+ 0+)]- 1 = G { £g A A \M 2u \\ (1) 

where G^J is the lepton phase-space integral, g A is the axial vector coupling 
constant and 



M 2v = (m e c 2 ) 



The Gamow- Teller (GT) matrix elements M 2 J are calculated by 



(2) 



m gt = v (o+llt+aiii^q+iif+ciiioj-) 

2v ^ \{Qpp + 2m e ci)+E N -E I ' [) 

where t is the isospin increasing/decreasing operator, er = 2s is the Pauli spin 
matrix, while Qpp is the Q value of the double-/? decay, and Ej and En are the 
energies of the initial and the intermediate states, respectively. The coefficient 
G^J is the lepton phase-space integral. Its values are given in Refs. [5J and [TD]. 
The Fermi (F) matrix elements M 2v are calculated by 

2 " n \(QPP + 2m e c2 ) + En - Ei ' 

The inverse half-life of the decay 0^ — > 2 J 

[rf; 2 (0+ -+ 2+)]- 1 = G^ 2 V|M 2 ,| 2 (5) 

is calculated in a similar fashion by M% v = (TO e c 2 ) 3 M^, T ' 2+ with 

Af GT,2+_ t /l V (2+||t+a||l+)(l+||t+a||0+) 
2v \l 3^(±(Qw+2m eC 2) + E N -E i y [> 

In this case, there is no Fermi contribution. The aim of this paper is the calcu- 
lation of the matrix elements M 2 J- , M 2v and Af^ T,2+ . 
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2.2 IBFM calculation of (3 decays 

The ingredients in the calculation are the matrix elements from even-even to 
odd-odd nuclei, (l^||t + cr||0f ), (0^ 1 1£ + 1 10^) , and from odd-odd to even-even, 
(0j||t + cr||l^), (0^||t + ||0^), which we now proceed to evaluate. 

A formulation of /3-decay in the proton-neutron interacting boson-fermion 
model (IBFM-2) was given years ago [4, Chap. 7] and [11] [12]. The microscopic 
theory gives the images of the Fermi and Gamow- Teller transition operators as 

o p = ^-V2JTT[p«p«] (0) , (7) 

3 

j'j 

where 

^ = -^'f;i'H|z§;i> 

The operator Pp stands for the boson-fermion image of the particle transfer 
operator. For the transitions from an even-even nucleus to an odd-odd nucleus, 
it can be either of the two operators: 

4^ = Qa] m + £ Gi' * f [d x a],}<£\ (10) 

y 

f 

where aj m is the fermion creation operator, s^ is the s-boson creation operator, 
and the z-component of d is related to the (i-boson annihilation operator by 
= ( — l) /i d_ jU .. In these operators, the distinction between the proton (ir) 
and the neutron (y) will be made later when necessary. The operator from an 
odd-odd nucleus to an even-even nucleus is 

M } = C*a jm + x MH\ (12) 

f 

4 W = OjfZim + E 9 n' $ x ^'1™ ' ( 13 ) 
i' 

where dj m is related to the fermion annihilation operator a,j m by dj m = (— l) 3- ™^- 
s is the s-boson annihilation operator, and d) is the d-boson creation operator. 



0) 
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The coefficients of the transfer operators are [1] 



Co 










Cm' 






Vj 1 











10 \ 1/2 1 



N(2j + iOf; 



10 ^ /2 



2j + lJ KK'> 



(14) 

(15) 
(16) 

(17) 



where u 3 and Vj are BCS unoccupation and occupation amplitudes, and the 
quantities K 7 K'j, K", 

/ \ 1/2 

. (18) 

if 

, / o ^j^ 2 (even;0 ] f |(n ;i + l)n d |even;0+)S/^ 2 'j V /2 , . 

^ = ( v 1 + 2 l^J 7V(2iTl) ' (19) 



is" ( (even;0J"|n s |even;0j") 
*j = I 



\ 2 (cvcn;0+|n d |even;0+) E/ %j \ 1/2 ,„„x 
V^J 2JT1 ) ' 



are calculated from the expectation values of the s-boson and d-boson numbers, 
h a , hd, and 

Pj'i = + v r u j) Qfj ( 21 ) 

Qi'i = (l'y , \\Y^\\qj). (22) 

If the odd fermion is a hole, then Uj and «j are interchanged, and the sign of 
Pjij is reversed in Eqs. ([H) to ([20]) . 



3 Calculation for i28,i30 Te ^i28,i30 Xe 

3.1 Energy levels and electromagnetic properties 

In order to calculate 2vf3j3 decay, we use the wave functions of the initial, in- 
termediate and final nuclei, in the present case the wave functions of 128 ' 130 Te, 
I28,i30j an( j I28,i30x e They are obtained from a calculation of the energy levels 
and electromagnetic properties. For purposes of checking the accuracy of our 
approach, we also calculate energy levels and electromagnetic properties of the 
odd-even nuclei 129 I, 127 Te, 131 I, 129 Te. The entire set of nuclei we calculate is 
shown in Table Q] 
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Table 1: Nuclei related to the double-/? decay from 128 ' 130 Te. 



initial 


final 


intermediate (odd-odd) 


related odd-even 


52 ie 76 


54 Ae74 


'53175 =^|Te 76 + p-n 


^I=^«Te+p 

127 Te= 128 Te _ n 


52 ie 78 


54 Ae 76 


1 5 3 3 U l77 = 1 5 3 2 U Te 78 +p-n 


131l = 130 Te+p 
i2|Te= 130 Te-n 



Table 2: The neutron and proton boson numbers, N v , N„, and the IBM-2 
parameters for 128 > 130r T e and 128 ' 130 Xe. The values are from Ref. [13 . The 
parameters that are not given in this table are set to zero. 



nucleus 


N v 


N v 


(McV) 


K 

(MoV) 






6,?2 

(MoV) 


«3 
(McV) 


4 U) 

(MoV) 


(McV) 


i'2H Te 


3 


1 


0.93 


-0.17 


0.50 


-1.20 


0.24 


-0.18 


0.30 


0.22 


128 Xe 


4 


2 


0.70 


-0.17 


0.33 


-0.80 


0.24 


-0.18 


0.30 


0.00 


13U Tc 


2 


1 


1.05 


-0.20 


0.90 


-1.20 


0.24 


-0.18 


0.30 


0.22 


130 Xo 


3 


2 


0.76 


-0.19 


0.50 


-0.80 


0.24 


-0.18 


0.30 


0.22 



3.1.1 12 8.i30 Te and 128130 Xe in IBM-2 
The Hamiltonian in IBM-2 is 

H B =e d (n dv + n d „) + k {Q* ■ Qf) 

+ 2^2 ((44 - 44) • {dvSn - dvSv)) 

+ E & H >< d t] (K) ■ K X 4] W ) 
K=l,3 

+ \ E cl K \[dixdt}W.[d v xl)W), (23) 

K =0,2, 4 

where 

- 4^ + 4d„ + X v [dl x d„]< 2 >, (24) 
q* = 4^ + 44 + X7T [4 x 4] (2) - (25) 

We adopt the parameters given in 13] , as shown in Table [2] Tables [3] and [4] 
show some of the calculated energy levels and their comparison with data. The 
agreement is very good. The same conclusion applies to the electromagnetic 
transition rates and moments, not shown here for conciseness. 
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Table 3: Energy levels in 128 Te and 128 Xe by the parameters in [13] . 





128 Xe 






128 Xe 




spin 


exp. 


cal. 


spin 


exp. 


cal. 




(MeV) 


(MeV) 




(MeV) 


(MeV) 


of 


0.000 


0.000 


ot 


0.000 


0.000 


n 


0.743 


0.739 


n 


0.443 


0.421 


4 


1.497 


1.672 


4 


0.969 


1.017 


4 


1.520 


1.528 


4 


1.033 


1.025 


0+ 


1.979 


1.971 


4 


1.430 


1.614 


4 


2.164 


1.940 


0+ 


1.583 


1.482 








4 


1.604 


1.751 



Table 4: Energy levels in 130 Te and 130 Xe by the parameters in [13] . 





130 Te 










spin 


exp. 


cal. 


spin 


exp. 


cal. 




(MeV) 


(MeV) 




(MeV) 


(MeV) 


4 


0.000 


0.000 


ot 


0.000 


0.000 


4 


0.839 


0.877 


4 


0.536 


0.511 


4 


1.588 


1.565 


4 


1.122 


1.223 


4 


1.633 


2.017 


4 


1.205 


1.228 


n 


1.886 


2.233 


4 


1.633 


1.795 


4 


1.965 


2.337 


4 


1.794 


1.677 


4 


1.982 


2.598 


4 




1.959 


4 


2.139 


2.010 


4 


1.808 


2.063 
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3.1.2 129 131 I and i^^Te in IBFM-2 

The Hamiltonian for odd-even nuclei (IBFM-2) is given by 

H = H B + H F p +V p BF . (26) 

The boson Hamiltonian H B is the core Hamiltonian ( 128 Te and 130 Te in the 
present case). The symbol p refers to it (proton) or v (neutron) depending on 
the odd fermion. The fermion single-particle Hamiltonian is 

where Ej is the quasi-particle energy of the odd particle, while hj is the number 
operator. We adopt the single-particle energies in [HJ[I2], shown in Table [SJ 
The quasi-particle energies Ej p are calculated in the usual BCS approximation 
with gap A = Ylj\J~A. In this BCS calculation, we include both positive and 
negative-parity orbits. The terms V p BF are the interaction between the bosons 
and the odd particle: 

rBF 

' p 



|(2) --[4<x<v] (2: 



' 7 

+ Hcrmitian conjugate 
+ A^2h l h dpl . (28) 

i 

The symbol p' indicates the other kind of nucleon; e.g., p' = v when p = n. For 
the orbital dependence of the interaction strengths, we adopt the parametriza- 
tion of Refs. (tjIH]: 

Tij = (uiUj - ViVj) Qij T, (29) 



4, - -m* 



1/2 

A, (30) 



fp&jk + 1) 
where 

= (uiVj + ViUj) Qij (31) 



. . - (l . I f ||y( 2 )|U. I j . ) + 5 (,! 20 ,,.n 



(32) 



The definitions of the parameters A and T are the same as that in [TT]. The 
exchange interaction in ([28)) with the same A corresponds to the total d-boson 
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Table 5: Single-particle energies for Z = 53 and N — 75 taken from Ref. [TTl IT2"] . 



orbit 


O57/2 


1^5/2 


1^3/2 


2S1/2 


0hn/2 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


proton 


0.00 


0.40 


3.00 


3.35 


1.50 


neutron 


0.00 


0.60 


2.50 


2.10 


2.00 



Table 6: Boson- fermion interaction parameters by Dellagiacoma [TTj . 



parameter 


r 


A 


A 




(MeV) 


(MeV) 


(MeV) 


proton in 12y . ial I 


0.60 


0.20 


-0.30 


neutron in 127 ' 129 Te 


0.30 


0.10 


-0.40 



number conserving part of that in [TT] . The factors Uj and Vj are interchanged 
if the odd nucleons are holes. 

The calculation splits into positive and negative parity levels. The positive 
parity orbitals for the odd proton and odd-neutron are: 

0ff7/2, 1^5/2) 14/2, and 2si/ 2 (33) 

while the negative parity orbital is /in/2- For the calculation of positive parity 
levels reported here, we use the parameters of [IT], shown in Table [6] Table [7] 
shows some of the low-lying energy levels. The agreement between calculated 
and experimental levels is good. 

3.1.3 128 > 130 I in IBFFM-2 

The intermediate states in 128 . 130 I are described by the proton-neutron interact- 
ing boson- fermion- fermion model (IBFFM-2) [15]. Because the nearest closed 
shell has Z = 50 and N = 82, the nuclei 128 ' 130 I are described as a system of an 
IBM-2 boson core and an odd proton and an odd neutron: 

128,130 T 128,130 . / OA \ 

53 : 75,77 = 52 Te 76 ,78 + P - H. (34) 

The odd neutron is treated as a hole. The related nuclei are summarized in 
Table [T] We include the same orbitals as in the previous subsection with s. p. 
e. given in Table [5j The Hamiltonian is 

H = H*+H* + V™ + Hi + V™ + Vres. (35) 

The boson and the fermion Hamiltonian parameters are those given in the pre- 
vious sections. The last term is the residual interaction between the odd proton 
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Table 7: Energy levels in 129 I, 127 Te, 131 I and 129 Te. The parameters of the 
boson Hamiltonian are taken from Ref. [13] as shown in Table[2] The parameters 
in the boson- fermion interaction are from as shown in Table [5] 





129 J 






i27 Tc 




spin 


exp. 


cal. 


spin 


exp. 


cal. 




(MeV) 


(MeV) 




(MeV) 


(MeV) 


7/2+ 


0.000 


0.000 


3/2+ 


0.000 


0.000 


5/2+ 


0.028 


0.155 


1/2+ 


0.061 


0.054 


3/2+ 


0.278 


0.505 


5/2+ 


0.473 


0.464 


5/2+ 


0.487 


0.645 


3/2+ 


0.502 


0.468 


1/2+ 


0.560 


0.641 


1/2| 


0.623 


0.516 








7/2+ 


0.685 


0.508 








3/2+ 


0.763 


0.558 








5/2+ 


0.783 


0.554 




131j 






m Te 




spin 


exp. 


cal. 


spin 


exp. 


cal. 


(MeV) 


(MeV) 




(MeV) 


(MeV) 




7/2+ 


0.000 


0.000 


3/2+ 


0.000 


0.000 


5/2+ 


0.150 


0.200 


1/2+ 


0.180 


0.148 


3/2+ 




0.703 


5/2+ 




0.609 


5/2+ 




0.711 


3/2+ 




0.621 


9/2+ 




0.794 


7/2+ 




0.660 
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Table 8: Proton-neutron residual interaction in 128 > 130 I. 



v s 


Vaa 


Va-a-S 


V T 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


-0.1 








0.05 



and the odd neutron given as |15j 

Vres = A-KV s 5{r n -r v )8{r v - RQ)5{r v - R Q ) 

— (oV • <X„) 

+ 4tt V aaS (er )S(r 7I - r v ) 6(r„ - R ) S(r, 

+ V T 3 ~2 ( w ' ^) ■ 

\ 7TI/ / 

The matrix elements between two quasi-particles are connected to those between 
two particles as 

quasi-particlc 

= ( U j' u U j' n U ju U ju + v 3l V K V jv V jv)(j'ufir : > J\VREs\jvj*, J) 

xJ2(2J' + l)\ j »' f J ' \{jX;J'\V RES \j u j' n ;J'). (37) 

~J7 L. Jv Jtt' o ) 

The strengths of the delta interaction (Vs), the spin-spin interaction (V aa ), the 
spin-spin-delta interaction (V aa s) and the tensor interaction (Vt) are determined 
from a fit to the experimental levels. The adopted values are shown in Table |SJ 

By diagonalizing the Hamiltonian (|35[) we obtain the energy levels and wave 
functions. The energy levels are compared with the experimental data in Table 
[9]and Fig. [TJ The agreement is fair. The ground state spins 1 + in 128 I and 5 + 
in 130 I are calculated correctly. The only disagreement is in the location of the 
state 2f which is calculated too high. However, for the purpose of the present 
paper, only the location of 1 + states is important. The location of + states 
is also of some importance, but these appear at higher excitation energy (805 
keV) and therefore are not shown in Table [9] and Figure Q] 

The electromagnetic transition operators in IBFFM-2 are 

- ]T ^Y,( u r^- v r^)(f\\^y Yi2) \\jH x ^] (2) ( 38 ) 

/»=*■, f j'j 



-Rq) 

(36) 



10 



Table 9: Energy levels in 128 I and 130 I. 











130 1 




spin 


exp. 


cal. 


spin 


exp. 


cal. 




(MeV) 


(MeV) 




(MeV) 


(MeV) 


It 


0.000 


0.000 




0.000 


0.000 




0.027 


0.193 




0.040 


0.195 




0.085 


0.045 


it 


0.043 


0.012 




0.128 


0.141 




0.044 


0.015 


3 2 + 


0.152 


0.090 


3+ 




0.091 






0.041 




0.048 


0.097 


4 




0.254 


4 




0.293 


n 




0.278 


n 




0.325 



0.3 -(2,3,4)+- 



> 



0.2 



0.1 



0.0 



128i 



4+ . 

1,2,3+" 

(3)+- 
1,2+. 
3+ - 



/ 

4+ *q 



\ 



3+ / ^ 

K 

2+ ' 



1+ 



2+ 
4+ 

2+ 
3+ 

4+ 

3+ 

3+ 
5+ 

1+ 



0.3- 



> 



0.2 



0.1 



0.0 



130t 



1,2,3+= 
3+ 

3+ - 

2,3+ ■ 



4,5+ 
3+ 



4± 
2+ 



/ 
/ 
/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ /- 



2+ 
4+ 



3+ 
2+ 



4+ 
3+ 



5+ 



exp IBFFM 



exp IBFFM 



Figure 1: Energy levels in 128 > 130 l. The experimental data are from [TBI [T9 
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and 



T(M1) = ( 5 « + « 

~~ 7J§X^''^ +v i""j){3'\\gi,pl + 9s,ps\\j)[a], x a^A, 

p=K,V V j'j ' 

(39) 

where is the boson angular momentum, I is the fermion orbital angular 
momentum, and s is the fermion spin. The effective charges and other co- 
efficients are taken from |12j . with which the electromagnetic properties in 
odd- A nuclei are explained very well: = e] 3 = 0.12 eb, e c tf :7T — 0.405 e, 
e e s <v = 0.135 e, g® = 1.3/x N) g® = -0.2/z N , gz i7r = l/z N , &, i7r = 3.910^ N , 
gi,v = 0/iN) .9s. — — 2.678/^n, where the spin g-factors have been quenched by a 
factor of 0.7. With these operators, we can calculate electromagnetic transitions 
and moments in 128 I and 130 I. They are given in Tables [TU] and [TTJ Unfortu- 
nately, not much experimental information is available, with the exception of the 
magnetic moment of the ground state of 130 I. Our calculated value /i 5 + = 3.12 
is in reasonable agreement with the experimental value, fi 5 + = 3.349 (7). The 
extensive tables are shown as a reference for additional experiments, if feasible. 

Table 10: Electromagnetic moments, n, Q, and transitions B(M1), 
B(E2), in 128 I. Since there are many levels with unclear spin, it 
is difficult to determine the level order. Therefore, the calculated 
decay rates are always written in the decreasing order of spin. If 
the actual decays take place in the opposite direction, then the 
spin factors of the reduced transition matrix elements need to be 
adjusted. 



level (s) 




exp 


cal 






(Mn) 


(Mn) 


1*0-1) 






1.74 


t*(#) 






1.95 


A*(2£) 






3.06 


Ml) 






2.74 


Mt) 






2.47 


ti*l) 






2.36 


MO 






2.01 


Ml) 






3.02 


B(M1;2+ -> 
B(M1;2+ -> 
B(M1;2+ -> 
B(M1;3+ -> 
B(M1;3+ -> 


4) 
ll) 
21) 
21) 
2 2 + ) 




K) 
0.130 
0.107 
0.00157 
0.000340 
0.00831 
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Table 10: Electromagnetic moments, /i, Q, and transitions B(M1), 
B(E2), in 128 I. Since there are many levels with unclear spin, it 
is difficult to determine the level order. Therefore, the calculated 
decay rates are always written in the decreasing order of spin. If 
the actual decays take place in the opposite direction, then the 
spin factors of the reduced transition matrix elements need to be 
adjusted. 



level(s) 




exp 


cal 


B(M1 


;3j" -> 


2f) 


> 0.020 


0.0510 


B(M1 


;3+ -> 


2J) 




0.000214 


B(M1 


;3j -> 


3f) 


> 0.003 


0.0531 


B(M1 


;4+^ 


3f) 




0.126 


b[mi 


;4+^ 


3j) 




0.00477 


B(M1 


;4-> 


3f) 




0.407 


b[m\ 




3f) 




0.274 


B(M1 


;4+^ 


4) 




0.213 


B(Ml 




4) 




0.00118 


b[m\ 


;5+^ 


4) 




0.000721 








(eb) 


(eb) 








-0.216 


QVt) 






-0.106 


Q(2+) 






-0.291 


Q$i) 






-0.460 


Q(3+) 






-0.310 


Q{±i) 






-0.372 








-0.525 






-0.590 








(e 2 b 2 ) 


(e 2 b 2 ) 


B(E2 


2+^ 


if) 




0.0206 


B(E2 


2+^ 


if) 




0.00552 


B(E2 


2+^ 


2f) 




0.00388 


B(E2 


3+^ 


If) 




0.000656 


B(E2 


3+^ 


2f) 




0.0108 


B(E2 


3+^ 


2J) 




0.00738 


B(E2 




If) 




0.0118 


B(E2 


3+^ 


2f) 




0.0192 


B(E2 


3+^ 


2f) 




0.00222 


B(E2 


3+^ 


3f) 




0.00362 


B(E2 


4+^ 


z i J 




0.00276 


B(E2 


4+^ 


2f) 




0.00110 


B(E2 


4+^ 


3f) 




0.00828 


B(E2 


4+^ 


3+1 

°2 J 




0.0117 


B(E2 




2+1 




0.000148 


B(E2 


4^ 


2+1 




0.00000796 
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Table 10: Electromagnetic moments, /i, Q, and transitions B(M1), 
B(E2), in 128 I. Since there are many levels with unclear spin, it 
is difficult to determine the level order. Therefore, the calculated 
decay rates are always written in the decreasing order of spin. If 
the actual decays take place in the opposite direction, then the 
spin factors of the reduced transition matrix elements need to be 
adjusted. 



Table 11: Electromagnetic moments, fj,, Q, and transitions B(M1), 
B{E2), in 130 I. 



levcl(s) 



cxp 



cal 



B(E2; 4+ ->• 3^) 
B{E2;4+ -> 3+) 
£?(£2;4+->4+) 
B(E2;5f -> 3f) 
B(E2;5+ -> 3t) 
B(E2;5f -> 4+) 
j?(£2;5+^4+) 



0.00109 
0.00415 
0.00179 
0.000124 
0.00501 
0.0146 
0.00504 



transition 



exp 



cal 



On) 



On) 
1.97 
1.96 
2.98 
2.30 
3.21 
2.58 
2.05 
3.12 



M(3f) 

M4) 



3.349 (7) 



OnT 



OnT 



B(M1 
B(M1 
B(M1 
B(M1 
B(M1 
B(M1 
B(M1 
B(M1 
B(M1 
B(M1 
B(M1 
B(M1 
B(M1 
B(M1 




0.115 
0.162 
0.0134 
0.00875 
0.00254 
0.0461 
0.0282 
0.0153 
0.0426 
0.0136 
0.0674 
0.650 
0.151 
0.00248 
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Table 11: Electromagnetic moments, fi, Q, and transitions B(M1), 
B(E2), in 130 I. 



transition 


exp 


cal 


B(Ml;5j-->4j) 




0.000461 






(eft) 


(e&) 






-0.121 


0(2f) 




-0.113 


Q(2+) 




-0.221 


Q(3f) 




-0.220 


Q(3+) 




-0.285 


Q(4f) 




-0.262 


Q(4) 




-0.378 


Q(5 + ) 




-0.377 








(e^) 


B(E2 


2+ -> 1+) 
z i ^ x i J 




0.0131 


B(E2 


2 + 1 + ) 




0.00496 


B(E2 


2+ -> 2+) 




0.00224 


B{E2 


3+ 1+) 




0.00199 


B(E2 


3+ 2+) 




0.0120 


B(E2 


3+ -> 2+) 




0.000507 


B[E2 


3 + 1 + ) 




0.00170 


B{E2 


3+ -> 2+) 




0.00190 


B(E2 


^2 ~~ ^ 2 2 ) 




0.00808 


B(E2 


3+ -> 3+) 




0.00146 


B(E2 


4+ 2+) 




0.00112 


B(E2 


4 + 2+) 




0.000512 


B(E2 


4 + 3+) 




0.00836 


B(E2 


4 + 3+) 




0.00286 


B(E2 


4 + 2+) 




0.0000952 


B(E2 


4 + -> 2+) 




0.000491 


B(E2 


4 + 3+) 




0.000126 


B(E2 


4 + 3+) 




0.00601 


B(E2 


4 + 4+) 




0.00117 


B(E2 


5+ 3+) 




0.00000478 


B(E2 


5+ -> 3+) 




0.00253 


B(E2 


5+ 4+) 




0.0104 


B(E2 


5 + 4+) 




0.00143 



3.2 Single-/? decay 

Single-^ decay matrix elements for 128,130^^.128,130! and ^s.isoj^ias.isoxg 
can be calculated using the formulas of Sect. 12.21 In the first process, 128 Te — ► 
128 I, the parent 128 Te and the daughter 128 I nuclei have the same boson numbers 
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(Nfr = lj-^Vi/ = 3). Therefore, the operators of (|10l) are applicable: 

Pp=A]^'\ PW=4«. (40) 

In the second process, 128 I — > 128 Xe, the even-even nucleus 128 Xe has (AT T = 
2,7V,, = 4) and the odd-odd nucleus 128 I has (N n =1,N„ = 3). Thus both for 
protons and neutrons, the transfer operator involves the creation of a fermion 
and the annihilation of a boson, and the operators of (I13|) are applicable: 

PP=BW'\ P<j»=BW). (41) 

This applies to the case 130 Te -> 130 I -> 130 Xe, too. 

We denote by M(F) and M(GT) the Fermi and Gamow- Teller matrix ele- 
ments, and by B(F) and _B(GT) their squares. From these we can calculate the 
log/t values for j3~ and /3 + /EC transitions using [TB] 

fil 63 

^ B(F)+ ^(GT) ^ + 1) <42) 

where is the angular momentum of the initial nucleus. The results, using the 
value gA = 1.269 from neutron decay [T7], are shown in Tables [T2HT51 column 
3. Some experimental information [T51 [T^] is available for 128 > 130 l — ^ I28,i30rp e 
decay, also shown in the tables, column 2. One can see that the magnitude of 
the calculated matrix elements is much larger than observed, resulting in a much 
shorter life-time. This is a well-known effect, due to quenching of the Gamow- 
Teller strengths in heavy nuclei. The quenched effective values of the axial 
vector coupling constant, <M,cff, can be obtained by a comparison between the 
calculated and experimental log ft values. We consider first the EC//3 + decay 
128 I (if) -> 128 Tc (Of). Using the experimental value [TS] log ft = 5.049 (7), 
we extract 

B(GT)[EC] = 0.102 (2). (43) 

The IBM value is 

B(GT) [IBM] = 1.676. (44) 

The ratio J B(GT)[IBM]/ J B(GT)[EC] = h 2 = 16.4, gives the hindrance factor 
h = 4.05. Large values of hindrance factor are expected in this region [5D]. The 
value of <7A,eff is then obtained from 



(fief) 2 B{GT) [IBM] = S(GT) [EC] ' (45) 

1.269 

9a,cS,ec/p+ = = 0-313. (46) 

A similar extraction can be done from the f}~ decay 128 I (if) — ► 128 Xe (Of ). 
The measured log ft = 6.061(5) [TS] gives 

B(GT)\p~] — 0.0100 (1). (47) 



yielding 
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Table 12: The log 10 ft values of electron capture from 128 I to 128 Te. The ex- 
perimental data are from |18j . 

transition exp cal quenched 

if Ot 5.049 (7) 3.836 5.15 (9) 



Table 13: The log 10 ft values of j3 decay from 128 I to 128 Xe. The experimental 
data are from [18]. 



transition 


exp 


cal 


quenched 


lt- 


+ ot 


6.061 (5) 


4.665 


5.98 (9) 


lt~ 


■» 


7.748 (24) 


5.262 


6.57 (9) 


lf- 


+ 0+ 


7.84 (6) 


5.712 


7.02 (9) 


lt~ 


^2+ 


6.495 (7) 


5.212 


6.52 (9) 




+ 2+ 


6.754 (9) 


6.446 


7.76 (9) 



The IBM value is 

£?(GT)[IBM] — 0.248 (48) 

with ratio B(GT)[IBM]/B(GT)[/3-] = h 2 = 24.8 and h = 4.98, from which we 
obtain 

1 9fi9 

9A.eB.p- = — = °- 255 - ( 49 ) 

This is consistent within 10% with <7A,eff,EC • We adopt in the following the 
average value 

<ffi = 0.28(3), (50) 
where we have estimated the error 5 in the determination of <75Pcff by 



^ _ / (&4,cff,EC//3+ ~ a) 2 + (9A,ctf,p- ~ 9? ^ 

with ~g the average value 0.284. If we use this value, we obtain the results shown 
in Tables HMH column 4. 

There are no available experimental data for the EC transition 130 I (1J") — > 
130 Te (0^) and for the f3~ transition 130 I (if) -> 130 Xc (Of), and therefore it is 
not possible to extract gA.eS for these decays. If we assume that su.eff for 130 I 
decay is the same as for 128 I decay, we can calculate the values given in Tables 
[T4"l and [TBI column 4. We see here that the quenched value 8.9 for the transition 
130 I (bf) — > 130 Te (4f ) is in fair agreement with the experimental value [19] 9.5 
(2). This transition is highly retarded both in theory and experiment. 

The main purpose of this paper is, however, the calculation of all 1^ and 0^ 
states in the intermediate odd-odd nuclei and of the matrix elements M(GT)(0f - 
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Table 14: The log 10 ft values of electron capture from 130 I to 130 Te. 



transition 


exp cal 


quenched 


1+^0+ 

1^0+ 

If 0+ 


3.626 
5.809 
5.664 


4.94 (9) 
7.12 (9) 
6.98 (9) 



Table 15: The log 10 ft values of f3 decay from 130 I to 130 Xe. The experimental 
data are from [H)] . 



transition 


exp 


cal 


quenched 


5f- 


+ 4+ 


9.5 (2) 


7.623 


8.94 (9) 


5f- 


+ 4+ 


8.2 (1) 


7.698 


9.01 (9) 




+ 4 


8.7 (1) 


9.117 


10.43 (9) 




+ at 




4.838 


6.15 (9) 




-> 0s" 




4.993 


6.31 (9) 


it- 


^0+ 




5.903 


7.22 (9) 



l]y) and M(F)(0f — > 0^). We have calculated GT and F matrix elements to 
all states up to an excitation of 3 MeV in 128 > 130 I. The GT matrix elements 
and energies £^(1^) up to 1.5 MeV excitation energy are given in Tables [T6l 
and 1171 Those for excitation energy in the entire range 0-3 MeV are shown in 
Figures [2] and [3] The F matrix elements and energies £^(0^) up to 1.5 MeV 
excitation energy are given in Tables [18] and 1191 Those for excitation energy in 
the entire range 0-3 MeV are shown in Figures 0] and \5\ There are 142 1 + states 
in 128 I, 108 1+ states in 130 I, 53 0+ states in 128 I, and 40 0+ states in 130 I, up 
to this energy. The main features of our calculation are: (1) The GT + matrix 
elements are distributed almost uniformly over the entire region 0-3 MeV, for 
i28430j ^ i28,i30 Te ^ (g) Qn the contrary, the GT" matrix elements 128 ' 130 I <-> 
128 ' 130 Xe are concentrated in only one state. (3) The F + matrix elements are 
concentrated in few states in the energy range 2-3 MeV. (4) The F~ matrix 
elements are uniformly small. 

The strength distribution £?(GT+) = |M(GT)+| 2 for 12 8,i30 Te _^ 128,130] has 
been measured recently by Puppe et al. [21] by means of the ( 3 He, t) reaction. 
The behavior of the experimental strength distribution is in agreement with the 
calculated behavior, i.e., the strength appears to be almost uniformly distributed 
as in our Figures [2] and [3l The authors of [21] extract also the values 



(52) 



3 I S(GT)£ s Ho '*) : 0.079 (8), = °- 829 ( 50 ) 
3 I S(GT)£ s Ho '*) : 0.072 (9), ^ = 0.746 (46) 
where represents the sum of the strength up to 3 MeV. The value 0.079 (8) is 



128] 
130 T 
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Table 16: The GT matrix elements ( 128 Xe| 128 I) and ( 128 I| 128 Tc) and the energies 
of the 1+ states (up to E x = 1.5 MeV) in 128 I. 



N 


( i28 Xe| i28 I) 
(0+||i+a||l+) 


( W8 Ip 8 Tc> 
(l+||*+a||0+) 

\ J V ' ' ' ' J. / 


128 T 
(MeV) 


1 


0.4981 


1.2946 


0.0000 


2 


0.0026 


0.0718 


0.4422 


3 


0.0821 


0.6681 


0.4953 


4 


0.0650 


-0.2570 


0.5434 


5 


-0.0136 


0.2809 


0.6226 


6 


-0.0643 


-0.6161 


0.8016 


7 


0.0156 


-0.2266 


0.8303 


8 


-0.0045 


-0.1336 


0.9111 


9 


0.0450 


0.4246 


0.9810 


10 


-0.0944 


-0.3538 


1.0062 


11 


-0.0761 


-0.8158 


1.0512 


12 


-0.0257 


0.1285 


1.0603 


13 


-0.0749 


-0.3665 


1.1175 


14 


-0.0429 


-0.2305 


1.1790 


15 


0.0138 


0.0108 


1.2021 


16 


-0.0155 


-0.1304 


1.2445 


17 


0.0043 


-0.0991 


1.2508 


18 


0.0134 


0.2683 


1.2796 


19 


0.0133 


0.2790 


1.2941 


20 


-0.0334 


-0.5122 


1.3383 


21 


0.0160 


0.4144 


1.3825 


22 


-0.0580 


-0.6427 


1.4118 


23 


0.0445 


0.2464 


1.4373 


24 


0.0141 


-0.0476 


1.4696 
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Table 17: The GT matrix elements ( 130 Xe| 130 I) and ( 130 I| 130 Te) and the energies 
of the 1+ states (up to E x = 1.5 MeV) in 130 I. 



N 


( Xel I) 
(Qt\\t + (j\\lt) 

\"1 1 1" w 1 1 AT/ 


/13UJI131W 
(ltllt+crllOT") 


130 T 

EJlt) 

(MeV) 


1 


0.4085 


1.6473 


0.0125 


2 


0.0094 


0.2874 


0.6142 


3 


0.0096 


0.266 


0.6699 


4 


-0.0153 


-0.1761 


0.6962 


5 


0.0024 


0.4817 


0.8129 


6 


-0.062 


-0.2209 


0.8893 


7 


0.0156 


0.3474 


0.9711 


8 


-0.045 


-0.2679 


1.0236 


9 


0.0698 


0.8778 


1.0726 


10 


0.0166 


-0.4301 


1.1387 


11 


0.0173 


-0.0084 


1.171 


12 


-0.0326 


-0.4325 


1.209 


13 


0.0141 


0.0357 


1.2359 


14 


-0.0214 


-0.0663 


1.2544 


15 


0.0063 


-0.3612 


1.2993 


16 


0.0056 


0.309 


1.3876 


17 


-0.0096 


0.0135 


1.4949 



Table 18: The F matrix elements ( 128 Xe| 128 I) and ( 128 I| 128 Te) and the energies 
of the 0+ states (up to E x = 1.5 MeV) in 128 I. 





( i28 Xc| i28 I) 


( i28 I| 128 Te) 




N 


<0+|t+|0+> 


(o^|t + |o+> 










(MeV) 


1 


-0.0017 


0.0836 


0.6673 


2 


-0.0497 


-0.1303 


0.7731 


3 


0.0191 


0.1871 


0.9547 


4 


0.0044 


0.1475 


1.1158 


5 


-0.0256 


0.1093 


1.3182 


6 


-0.0174 


0.3785 


1.3382 


7 


-0.0319 


0.6042 


1.4194 


8 


0.0041 


-0.0838 


1.4338 
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Figure 2: The values of (l^||i+oj|0^) (top), (0f||<+cr||l^) (center) and 
(0+||t+ ( 7||l+)(l+|K+a||0+)/(i(Q / 3 /3 + 2m e c 2 ) + E N - Ej) (bottom), for the 
double-/? decay from the lowest + in 128 Te to the lowest + in 128 Xe through 
the intermediate 1 + in 128 I, plotted as a function of the excitation energy of 1 + . 
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Figure 3: Same plots as Fig. [2] but for 130 Te decay. 
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Figure 4: The values of (C)£||i + ||0|) (top), (0|||t+||0^) (center) and 
(0f\\t+\\0+)(0+\\t + \\0i)/(±(Q l 3p+2m e c 2 )+E N ~E I ) (bottom), for the double-/3 
decay from the lowest + in 128 Te to the lowest + in 128 Xe through the inter- 
mediate + in 128 I, plotted as a function of the excitation energy of + . Note 
the small scale in the bottom panel of the figure. 
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Figure 5: Same as Fig. 0] but for 130 Te decay. 
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Table 19: The F matrix elements ( 130 Xe| 130 I) and ( 130 I| 130 Te) and the energies 
of the 0^ states (up to E x = 1.5 MeV) in 130 I. 





( 130 Xe| i3u I) 


(i30 T po Xe) 




N 


<0+|i+|0+> 


<0+|* + |0+) 










(MeV) 


1 


0.0245 


0.0756 


0.8052 


2 


-0.0095 


0.0589 


0.9045 


3 


-0.0477 


0.0466 


1.1851 


4 


-0.0136 


0.2426 


1.2933 


5 


0.0002 


-0.0670 


1.4032 



in fair agreement with 0.102 (2) obtained from dSHH] and the value 0.087 (10) 
calculated by the authors from the data of [22]. From a comparison between 
the values in ([52]) and the calculated IBM values 

128 I fl(GT) : 1.676, V = 15.09 

(53) 

130 I B(GT) : 2.714, ^ = 16.22 

we can extract the values of SU,eff,( 3 He,t)- We obtain, for 128 Te (Of ) — !• 128 I (1^), 
gA eS = 0.276 and for 130 Te (0+) 130 I (1+), g A , eS = 0.207. The values so 
extracted are smaller than those extracted from EC//? + and /3~ decay, especially 
for 130 Te. This may have to do with the way in which J B(GT)( 3 He, t) is extracted 
from the cross section, causing a tension between i?(GT)[EC] = 0.102 (2) and 
B(GT)[ 3 He, f] = 0.079 (8) for 128 I (1+) o 128 Te (0+). 

If we use the value <ju,eff = 0.28 (3) obtained from EC//3 + and j3~ decays, 
we calculate the quenched values 

128 I B(GT) [IBM-quenched] : 0.082, V [IBM-quenched] = 0.735 

" (54) 
130 I B(GT) [IBM-quenched] : 0.132, ^[IBM-quenched] = 0.790. 

The values for 128 I are in good agreement with the ( 3 He, t) values, but only 
in fair agreement for 130 I, where we overestimate the summed strength by 5%. 
We conclude that IBM-quenched calculated values with a single value gA,eS — 
0.28 (3) provide a good description of all experimental data EC//3 + , /3~ and 
( 3 He, t) in 128 I decay, and a fair description of ( 3 He, t) in 130 I decay. 

There is no reported measurement of the strength distribution B(F + ) = 
|Af(F)+| 2 . In Ref. [21] only the IAS 0+ state at 11.948 MeV in 128 I and at 
12.718 MeV in 130 I is identified. However, some excess strength in the region 2- 
3 MeV, especially around E x ~ 2.7 MeV where we predict the F + strength to be 
concentrated, is seen in Fig. 1 of [21] . It would be of great interest to investigate 
this point further, since it will clarify the question of isospin violation in this 
mass region. 
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Table 20: Nuclear matrix elements M^ T , Mj„ of transitions from the ground 
state of 128 < 130 Te to some states in 128430 Xe. The sign of M^J is chosen to be 
positive. 



transition 


12« Te ^ 12 » Xe 


l30 Te ^l30 Xe 


GT 






0+^0+ 


0.297 


0.273 


0+^2+ 


0.00718 


0.00639 


0+^0+ 




0.668 


F 






0+^0+ 


-0.0353 


-0.0309 


0+ 0+ 




-0.112 



3.3 Double-/? decay 

The individual matrix elements of t + a and t + are then combined as in Eqs. ([3]), 
((4]) and (J6j> to obtain the matrix elements for 2^/3/3 decay. In evaluating the 
denominators in Eqs. ([3]), Q and ©, the following experimental Q- values are 
used: Qp (gs) ( 128 Te -> 128 I) = -1.2518, Q w (gs) ( 128 Te -> 128 Xe) = 0.868, 
Q/3 (gs) ( 130 Te -J- 130 I) = -0.420, Q (ifi (gs) ( 130 Te -J- 130 Xe) = 2.529, in units of 
MeV. When calculating decay to 2+ states, we use the experimental values Qpp 
(2+) ( 128 Te -> 128 Xe) = 0.425, Q pf3 (2+) ( 130 Te -> 130 Xe) = 1.993, also in MeV. 
The calculated values of En (exc) in the intermediate odd-odd nucleus and of 
the individual matrix elements are then combined and summed to give the value 
of the nuclear matrix elements shown in Table [20l The Fermi matrix elements 
in this table are a factor of approximately 10 smaller than the Gamow- Teller 
matrix elements. Introducing the quantity xf = M^/M^J, of Ref. [9], we find, 
for 0+ 0+ transitions xf = -0.119 in 128 Te decay and xf = -0.113 in 130 Te 
decay. If isospin was a good quantum number, xf = 0. Table [50] indicates 
that there is a small isospin violation in our wave functions. The individual 
contributions to the sums are also shown in the bottom panels of Figures [2HS1 
for i28,i30 Te ^i28430 I ^i28430 Xe _ We find that the GT sum jg dominated by the 

lowest if state in 128 > 130 l. This is the most important result of this paper. Our 
calculation is consistent with (1) the single state dominance (SSD) hypothesis 
[H [Ml 123 and (2) the Fermi-surface quasi-particle model of Ejiri [H] . The F 
sum is very small and receives most of its contribution from two states at E x = 
2.179 and 2.701 MeV in 130 I, and at E x = 2.282 and 2.443 MeV in 128 I. 
The matrix elements 



\M 2 A=g\ 



(55) 



for 0^ — > Of decay can be extracted from experiment using the observed half- 

r l/2 



life t?J 2 , Eq. dU). The extracted values are [TU] 0.044 (6) for 128 Te decay and 



0.031 (4) for 130 Te decay. These values should be compared with the calculated 
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Table 21: Two-neutrino double-/3 decay matrix elements, \Mi V \ in IBFM. 



exp calc quenched 
Te 0.044 (6) (1514 0.040 (8) 
Te 0.031 (4) 0.470 0.037 (8) 



values 



M? C \=9 2 A 



M GT_f9vYMF 
V OA ) 



L 2u \ I al 2v 



(56) 



From the values in Table |2D] and g A = 1.269, g v = 1, we have 0.514 for 128 Te 
and 0.470 for 130 Te decay. Under the assumption that gA is quenched to gA,cS,/3/3 
while (gv/gA) is not, thus writing 

|Mr nched |=5i, eff ,^l^ alc | ) (57) 

we can extract gA,eft,pp = 0.293 for 128 Te and gA,eff,0/3 = 0.257 for 130 Te, in fair 
agreement with the values extracted from single-/? decay, gAeS p = 0.313 for 

128 Te ^12 8l and o 255 for 12 8l ^128 X ^ 

Assuming that the quenching of gA is the same in both single-/3 and 2vf3f3 
decay and using the adopted value gA, c s = 0.28 (3) of Eq. (|50|) we obtain 
the quenched values in Table [2P1 in reasonable agreement with experiment. It 
appears from this table that knowledge of single-/? decays allows one to reliably 
calculate 2v/3f3, as emphasized by Ejiri [55], and thus predict the 2^/3/3 half-life 
in cases where it has not been measured. This statement, however, relies on 
our assumption leading to Eq. (|57l) . Quenching of matrix elements in a given 
model calculation arises from two effects: (i) The limited model space in which 
the calculation is done and (ii) coupling to non-nucleonic degrees of freedom 
(A, N*, ...). For the second part we expect gA to be quenched and gv not, 
due to the conserved vector current hypothesis (CVC). For the first part, it is 
reasonable to expect that both gA and gv be quenched. While for gA there are 
experimental data to extract gA,cS from single-/? decay, as we have done in Sect. 
13. 2\ there are no data to extract gv,oS, and thus our assumption that gv/gA in 
(|57[) is unquenched is speculative. 



(0) 

2u,0+ ->0 



Using the values in Table [2H and the phase space factor of [TU] , G 
(i28 Te ^128 Xe) = o 26Q x 10 -2i yr -i and G^° ) 0+ _ f0+ ( 130 Te -^ 130 Xe) = 

1529 x 10- 21 yr" 1 , we calculate from r" 1 = G { °J |M 2 q " cnchod | 2 , the half-lives 
r( 128 Te) = 0.23 (9) x 10 25 yr and r( 130 Te) = 0.48 (19) x 10 25 yr to be compared 
with the experimental values 0.19 x 10 25 yr ( 128 Te) and 0.68 x 10 21 yr ( 130 Te). 

In addition to matrix elements to 0^~ , we have also calculated matrix elements 
to 0+. This state is located at 1.583 MeV in 128 Xe and at 1.793 MeV in 130 Xe. 
Therefore Q w (0+) ( 128 Te -> 128 Xe) = -0.715, Q w (0+) ( 130 Te -> 130 Xe) = 
+0.736. Decay to 0^" is possible for 130 Te decay, while for 128 Te decay it is not. 
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The values of the GT and F matrix elements to 0^ are also shown in Table 
[20l They are larger than those to 0^~ due to the smaller energy denominator in 
Eqs. ([3]) and (j3|. They can be combined as in Eq. (I55|) to give 1.187. Using the 
same quenching factor gA,cS — 0.28 (3) as before, we obtain M| !y ' qucnchcd (0 ] f — > 
0%) = 0.093 (18). 

2^/3/3 decay to 0^" has not been observed in 130 Te decay. It has so far been 
observed only in 100 Mo and 150 Nd decay Q]. Our calculation indicates that it 
may be observed. Using the values of Ref. [10], for + + = 1529 x 10~ 21 

yr" 1 and G^ 0+ ^ 0+ = 0.0757 x RU 21 yr" 1 , and the values |M 2 f ' quonc hod (0+ -> 
0+)| - 0.037 and |M 2 f' quenc hcd (0+ ->■ 0+)| = 0.093 we find 

T 2 ,( 13 °Te(0+)^ 130 Xe(0+)) , N 

\, n — — = 3200. 58 

T2u( 3 °Te(0+)^ 1 3o Xe(0 +)) 

The ratio is of the same order of magnitude of 

M^mof) MQf a _ 2790 . (69) 

n*{ Te(0f) ->«° Xe(0+)) 
It is however far larger than in 100 Mo decay where the observed ratio is p] [10] 
T 2i ,( 10 °Mo(0+) ^ 10 ° Ru(0+)) 



r 2i ,( luu Mo(0+)^ 130 Ru(0+)) 
and thus it may be difficult to observe. 



70 (10) (60) 



3.4 Sensitivity to parameter assumptions 

Single-/? and double- /3 decays are particularly sensitive to the occupation prob- 
abilities Uj, Vj of single-particle orbits, as one can see from Eqs. (| 14 1) - ([!?)) . To 
test this sensitivity, we have redone the calculation with another set of single- 
particle energies, as proposed by Fujita and Ikeda [20J, shown in Table |2"21 This 
set is rather different from that in Table [5] most notably by the location of the 
O57/2 level and by the inclusion of the 0/ig/ 2 proton and 0gg/ 2 neutron levels. 
It does not reproduce accurately spectra of odd-even and odd-odd nuclei in the 
region, but it is considered here to test the sensitivity to the choice of single 
particle energies. 

With this set we re-calculate the distribution of matrix elements for 128 i 130 Te— >1 
and 128,130 I— >Xe, as shown in Figures [BJ and [7] for GT. By comparing with Fig- 
ures and [3] one can see that the qualitative features, including the single-state 
dominance, are unchanged. The distribution in the leg 128 > 130 Te— >l is prac- 
tically unchanged. However, a major change occurs in the magnitude of the 
GT matrix elements from 1 + in 128 < 130 l to 128,130 Xe, which are smaller than 
in the calculation with the s. p. e. of Table [5] This change results in a reduc- 
tion of the double-/3 matrix elements |M2„| as shown in Table |2"51 The same 
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Table 22: Single-particle energies taken from Ref. 20 . 



orbit 


O.gg/2 


1^5/2 


0#7/2 


2si/ 2 


1^3/2 


0hn/2 


O/19/2 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


proton 




0.00 


-0.63 


2.40 


2.30 


2.20 


7.70 


neutron 


-4.08 


0.00 


0.42 


1.90 


2.20 


2.40 





Table 23: Nuclear matrix elements M21, of the transitions from the ground 
states of 128 > 130 Te to some states in 128 > 130 Xe calculated with the single-particle 
energies of [20] . 



transition 


12« Te ^12« Xe 


130 Te ^130 Xe 


GT 






0+^0+ 


0.178 


0.152 


0+^2+ 


0.00334 


0.00201 


0+->0+ 




0.557 


F 






0+^0+ 


-0.0297 


-0.0268 


0+ -> 0+ 




-0.110 



situation occurs for the F matrix elements. By repeating the same procedure 
as in Section [3.21 one can extract the values of gA,eg,p and <L4,eff,(3He,t) f° r this 
set of s. p. e.. For 128 I (l+)-> 128 Te (Of) decay, g A off = 0.292, and for 128 I 
(l+H 128 Xc (0+) decay, g A , e3 , p = 0.413. For 128 Te (0+)-> 128 I (if) ( 3 He, t), 
one has gA,eS,( 3 He,t) — 0.257. These values are somewhat inconsistent with each 
other. Taken on the average, they lead to a larger value of gA,eS — 0.321. By 
repeating the same procedure as in Section [3.31 one can also extract the values 
g A ^ M = 0.338 ( 128 Te) and 0.373 ( 130 Te). The extracted values of g A , c s,f3 are 
here also somewhat inconsistent with those extracted from single-/? decay and 
( 3 He, t). Nevertheless, assuming gA,eff,2v/3i3 = gA,oB,/3, and using the average 
value gA,eS = 0.35 (6) from single-/3, we obtain the values in Table IIMl These 
are in reasonable agreement with experiment, although the agreement is not as 
good as with the s. p. e. of |11) . 

3.5 Sensitivity to truncation to E x < 3 MeV 

The calculations reported in the previous subsections are based on contributions 
of states in the intermediate odd-odd nucleus with E x < 3 MeV. It is of interest 
to investigate how likely or unlikely is that states above 3 MeV could contribute 
significantly to two-neutrino decay. 

The 1 + and + states with E x > 3 MeV are built from two contributions: 
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Figure 6: The values of (l^||t+cr||0^) (top), (0^ \\t+ a\\l+) (center) and 
(Oi\\t+a\\l+)(l+\\t+a\\0+)/(^Q^ + 2m e c 2 ) + E N - Ej) (bottom) for the 
double-/3 decay from the lowest + in 128 Te to the lowest + in 128 Xe calcu- 
lated by the single-particle energies of [2D] , 
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Figure 7: The same plots as Fig. [3] for the decay 130 Te-> 130 Xe through 130 I 
calculated by the single-particle energies of • 
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Table 24: Two-neutrino double-/? decay matrix elements, \M2 U \ in IBFM with 
single particle energies of [20] . 



exp 


calc 


quenched 


128 Te 0.044 (6) 
130 Te 0.031 (4) 


0.315 
0.271 


0.039 (13) 
0.033 (11) 



(1) states constructed from single-particle orbitals within the model space of 
Table [5] (2) states constructed from single-particle orbitals above the shell gap 
at 82 or below the shell gap at 50. To investigate the contribution of states 
with E x > 3 MeV within the model space of Table [5j one can simply extend the 
calculation from the current lowest ~ 100 1 + states and ~ 50 + states to larger 
numbers. It appears that the properties of the strength distributions remain 
the same and that therefore states of this type will not contribute significantly 
to two-neutrino decay. States constructed from single-particle orbits outside 
the model space of Table [5] will appear in the spectrum at energies above the 
shell gaps, E x > 5 MeV. We have in fact investigated their contributions by 
including the proton orbit Oft.9/2 and the neutron orbit Ogg/2, as in Table l22l 
and concluded that also this type of states will not contribute significantly to 
two-neutrino decay. The argument is as follows. Inclusion of excitations across 
major shells will give major contributions to the strength distribution for the 
GT~ and F~ "legs", Te — ► I, especially in the region E x ^ 10 MeV, where 
the giant GT resonance (E x ~ 14 MeV) and the Isobaric Analogue State, IAS 
(E x ~ 12 MeV) are located. However, because of their composition in terms of 
single-particle states, we expect its contribution to the GT + (or F + ) "leg", I — > 
Xe, which will still remain concentrated in few low-lying states, Figs. [5] and 
center panel. It is therefore, in our opinion, quite unlikely that states above 3 
MeV will contribute significantly to two-neutrino decay. 



4 Conclusion 

In this article, a detailed investigation of the ten nuclei: 128 . 130 Te, 128 > 130 I, 
I28,i30x e ^ I29,i3ij anc [ 127,129^ w fthin the framework of the interacting boson 
model-2, IBM-2, and its generalizations IBFM-2 and IBFFM-2 has been done. 
The parameters needed in this investigation have been obtained as much as 
possible from the available experimental information. The wave functions so 
obtained have been used to calculate single-/? and 2v(5j3 matrix elements. 
The main results of our investigation are: 

(1) The mechanism of 2v/3f3 in these nuclei appears to be single-state dominance 
(SSD). 

(2) The GT strength in 12f M30 Te ( Q +) _> 128,130j ig even i y distributed. 
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(3) The GT strength in 128 < 130 I (1+) -> i28,i30 Xe ( Q +) is con centrated in one 
state. 

(4) Use of a single value gA. c s = 9A,eS,/3 = 9A,eS,0/3 appears to describe well 
both single-/? and 2f/3/3 decay. 

(5) The results are very sensitive to the choice of single-particle energies, most 
notably the weak branch 128 ^ 130 l _j.i28,i30 ^e. However, when the renor- 
malization of gA is taken into account through fitting the single-/3 decay, 
different choices of s. p. e. give similar results for 2f/3/3, but with varying 
degree of accuracy. The best choice appears to be that of the s. p. e. of 
[TT] which describes all observed quantities fairly: energies, electromagnetic 
transitions and moments, single-/? matrix elements, ( 3 He, t) strength distri- 
butions, and double-/? matrix elements, in the quenched approximation. 

Our best estimates of 2f /3/3 matrix elements are therefore those given in Table 
HU |Af 2 c 5| = 0.040 (8) for 128 Te(0f) ^ 12S Xe(O^) and \M^\ = 0.037 (8) for 
130 Te(0^) ^ 130 Xe(O^), and of the 2v0P half-lives r ( 128 Te) = 0.23 (9) x 10 25 
yr, t ( 130 Te) = 0.48 (19) x 10 21 yr. 

Our extracted values gA, e s = 0.28 (3) and su.eff = 0.35 (6) for the single- 
particle levels of Table 151 and l22| respectively, are rather low. The values of gA.cS 
depend on mass number, A, and on the nuclear model used in their extraction. 
A preliminary study of gA.cB m the Interacting Boson Model (IBM-2) in the 
closure approximation for 2f/3/3 decay and the Interacting Shell Model (ISM) 
has been done in [30j . It has been found that gA,eS has a smooth dependence that 
can be parametrized as gA.cff — 1.269A -018 plus shell effects. The extracted 
values of gA. off in the mass region, A ~ 130, are gA.cS — 0.5 for IBM-2 and 
~ 0.6 for ISM. A similar analysis has been done within the framework of QRPA 
[31] with similar results. We intend to continue the study of gA,eS within the 
framework described in the present article to understand how general is the 
result presented here, and also to study the related question of the extent to 
which gv is quenched in heavy nuclei, if at all. 

The question of the impact of the small value of gA, off found in /3-decay 
and 2f/3/3-decay to Of/3/3-decay is the subject of much debate. While only 
GT (1 + ) and F (0 + ) multipoles contribute to allowed /3- and 2f /3/3-decay, all 
multipolcs (l + ,2~,3 + , ...), (0 + ,l~,2 + , ...) contribute to Of/3/3 decay. It is not 
clear whether or not the higher multipoles are quenched. Nonetheless, since 
1 + and + still provide the largest contributions, we expect the quenching of 
gA and gv to play an important role in Of/3/3 decay. Since, in view of the fact 
that gA appears to the fourth power in the decay rate, the quenching of gA has 
major repercussions on Of/3/3 experiments, we plan to investigate this problem 
in depth in subsequent papers. 
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